Modulation of a quantized vortex street with a vibrating obstacle 
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Dynamics of a superfluid flow past an obstacle are investigated by solving the Gross-Pitaevskii 
equation numerically. For an appropriate velocity and size of the obstacle, quantized vortices are 
periodically generated in the wake, which form a Benard-von Karman vortex street. It is found 
that vibration of an obstacle modulates the vortex street breaking a symmetry. 



I. INTRODUCTION 

When a cylindrical obstacle is moved in a fluid with 
an appropriate velocity, vortices are periodically shed in 
the wake, forming a long-lived vortex street. This phe- 
nomenon was experimentally studied by Benard [J] and 
the stability of the vortex street was theoretically investi- 
gated by von Karman [2, 3] in the early twentieth century. 
The Benard-von Karman (BvK) vortex streets are seen 
in a wide scale throughout nature: from a bathtub to the 
atmosphere of a planet. 

The behaviors of a classical viscous fluid past an obsta- 
cle are characterized by the Reynolds number Re, which 
is proportional to the size and velocity of the obstacle 
and inversely proportional to the kinematic viscosity. A 
stationary flow is stable for a small value of Re. When Re 
exceeds a critical value, an instability sets in and vortices 
are shed in the wake, forming the BvK vortex street. For 
a very large value of Re, the flow behind the obstacle be- 
comes turbulence. Thus, in a classical fluid, the viscosity, 
which is included in Re, plays an important role for the 
vortex street generation. 

In contrast, Sasaki et al [4] theoretically predicted 
that a BvK-like vortex street can also be generated in 
superfluids, in which viscosity is absent and a vortex is 
quantized. This prediction has been made by numeri- 
cal simulations of the Gross-Pitaevskii (GP) equation, 
and therefore the phenomenon can be observed in the 
systems that are described by the GP equation, such 
as an atomic-gas Bose-Einstein condensate [B| and an 
exciton-polariton superfluid in a semiconductor micro- 
cavity [6|, l7| . Theoretical studies on the dynamics of su- 
perfluids behind a moving obstacle have been performed 
by many researchers [8I-KL4}. 

In the present paper, we extend the study in Ref. [4] 
to the case in which a moving obstacle oscillates in the 
transverse direction, i.e., the obstacle meanders. In clas- 
sical fluids, it has been found that such a vibrating cylin- 
der generates various types of vortex streets, in which the 
symmetry between two rows of vortices is broken [TB|, 03 • 
We will show that such symmetry breaking occurs also 
for the BvK vortex street in a superfluid. When the oscil- 
lation frequency is close to the vortex shedding frequency, 
the vortex street is modulated and becomes asymmetric. 

This paper is organized as follows. Section [II] for- 
mulates the problem. Section IIII Al reviews the case of 
straight motion of an obstacle [4] and Sec. IIII Bl shows 



the case of sinusoidal motion, which is the main result of 
the paper. Section HVl gives conclusions to the study. 



II. FORMULATION OF THE PROBLEM 



by 



We consider the two-dimensional GP equation given 
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(i) 



where ifj(x,y,t) is the macroscopic wave function, 
V(x, y, t) is an obstacle potential, and g is an interaction 
coefficient. The obstacle potential is assumed to have a 
Gaussian shape as 



V{r,t) = V eM-[r-r po t(t)} 2 /d 2 }, 



(2) 



where Vq is the peak intensity, r pot (t) is the position of 
the center, and d is the 1/e width. For an atomic Bose- 
Einstein condensate, such a potential can be produced 
by a blue-detuned Gaussian laser beam. In the present 
study, we consider sinusoidal motion of the potential as 



r po t(t) = —vtx + asm(ut)y, 



(3) 



where x and y are the unit vectors, v > is the velocity 
in the x direction, and a and uo are the amplitude and 
frequency of the oscillation. 

We numerically solve Eq. (pQ) using the pseudospec- 
tral method [TtJ. The initial state is the ground state of 
Eq. (pQ) with the initial position of the Gaussian poten- 
tial, which is obtained by the imaginary-time propaga- 
tion method. The norm of the wave function, J \ijj\ 2 dr, 
is chosen such that the density no far from the obsta- 
cle potential becomes a desired value. To break the ex- 
act numerical symmetry, small white noise is added to 
the initial state. We take a large enough space and the 
boundary condition does not affect the result. 

The characteristic length scale of the system is the 
healing length defined by £ = %/ '(mgno) 1 / 2 . The sound 
velocity for the uniform density no is v s = (gno/m) 1 / 2 . 
The characteristic time scale and energy are defined as 
r = h/(gno) and e = gno, respectively. For a typical 
experimental system of an atomic Bose-Einstein conden- 
sate, £ ~ 0.1 /im, v s ~ 1 mm/s, and r ~ 0.1 ms. For an 
exciton-polariton superfluid in a semiconductor micro- 
cavity, £ ~ 1 /im, v s ~ 10 6 m/s, and r ~ 1 ps. 
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FIG. 1: Density (IV^ 2 ) and phase (arg ip) profiles for the straight motion of the obstacle, a = 0. (a) v = 0A8v s , d — 1.26£, and 
Vb = lOOe. (b) v — 0.44v s , d — 2.24£, and Vb = lOOe. The arrows in (a) indicate directions of propagation of vortex dipoles. 
The arrows in (b) indicate directions of rotation of corotating vortex clusters. The field of view is 250£ x 95£. See the movie [HI 
showing the dynamics in (a) and (b). 



III. NUMERICAL RESULTS 
A. Straight motion of an obstacle 

I 

I 



The GP equation in Eq. (pp) has a scaling property. 
Normalizing the length, time, and density by £, r, and no, 
respectively, we can eliminate the interaction coefficient 
g from the equation. The independent parameters are 
therefore Vb, <i, v, a, and u. 



First, we consider the case of a = in Eq. (j3j), i.e., 
the potential moves straight in the — x direction. It has 
been known that pairs of quantized vortex and antivor- 
tex, which we call vortex dipoles, are released behind 
an obstacle above some critical velocity Increasing 
the velocity, the quantized vortices in the wake form 
two kinds of typical patterns [4|, as shown in Fig. [TJ In 
Fig. IHa), periodically released vortex dipoles incline al- 
ternately breaking the symmetry with respect to the x 
axis. Since vortex dipoles propagate in the directions as 
indicated by the arrows in Fig. DJa), the inclined vor- 
tex dipoles propagate obliquely, and the vortex pattern 
spreads in the downstream. 

Figure QJb) shows an interesting vortex pattern, which 
was first reported in Ref. [4]. The building block of the 
vortex pattern is not a vortex dipole but a pair of coro- 
tating vortices with the same circulation, which we call a 
"vortex cluster". The vortex clusters with opposite cir- 
culations are released alternately, which form two rows in 
the wake. This vortex shedding dynamics and resultant 
vortex pattern is similar to the BvK vortex street in clas- 
sical fluids. As shown by von Karman 0, Q, the vortex 
pattern as generated in Fig. DJb) is very stable and has 
a long lifetime. 

From the movies showing the dynamics in Fig. [TJ it 
seems that two different instabilities are responsible for 



the vortex pattern formations. The first one is the insta- 
bility that inclines the vortex dipoles in Fig. [I] (a); At the 
first stage, the vortex dipoles have the symmetry with re- 
spect to the x axis, and after some time the inclination 
of the vortex dipoles starts. The second one is the in- 
stability that causes a transition from the vortex dipole 
generation to the BvK vortex street generation; In the 
movie corresponding to Fig. [1] (b), we can see that the 
former changes to the latter dynamically. 

For the velocity v slower than the critical velocity of 
vortex generation, the flow pattern is a stationary lam- 
inar flow. For a very large velocity, the flow pattern 
becomes irregular or a C her enkov- like pattern [l9|, [20j . 
The periodic vortex shedding as shown in Fig. [1] is re- 
stricted between these two parameter regions. (The de- 
tailed parameter diagram is obtained in Ref. [4].) The 
vortex dipole generation as in Fig. QJa) tends to occur 
for a smaller obstacle, and the BvK-like vortex street 
tends to appear for a larger obstacle. For an interme- 
diate size of an obstacle, the former (latter) occurs for 
a slower (faster) velocity, where the vortex dipole gen- 
eration changes to the BvK-like vortex street with an 
increase in the velocity. 

B. Sinusoidal motion of an obstacle 



We next consider the case of a ^ in Eq. ((3]), i.e., the 



obstacle potential follows a sinusoidal trajectory. The 
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FIG. 2: Density (l^l 2 ) and phase (arg ip) profiles for the sinusoidal motion of the obstacle, (a) a = 0.32£ and u — O.lr -1 . (b) 
a — 0.64£ and u = O.lr -1 . Other parameters are the same as those in Fig. [T] (b). The field of view is 250£ x 95£. See the 
movie [18[ showing the dynamics in (a) and (b). 



vortex patterns generated by the obstacle are shown in 
Fig. [2j where the parameters are the same as those in 
Fig.QJb) except the value of a. The frequency uj is almost 
the same as the shedding frequency of vortex clusters. 
The vortex pattern in Fig. [2^a) is similar to the BvK- 
like vortex street in Fig. QJb). However, we find that 
the symmetry between the two rows of vortex clusters 
is broken, that is, the distance between two vortices in 
each vortex cluster in the upper row is larger than that 
in the lower row. The frequency of corotation is inversely 
proportional to the distance between the two vortices in 
each vortex cluster. The behaviors of the upper and lower 
rows can be exchanged in the dynamics, which is seen in 
the movie of Fig. EJa) [H|. 

The asymmetry is enhanced for a larger oscillation am- 
plitude a, as shown in Fig. Efb). The size of the vortex 
clusters in the upper row becomes very large as they flow 
to the downstream, making the vortex pattern different 
from the BvK vortex street. In the movie of Fig. Efb), 
it seems that the exchange between the upper and lower 
rows, that is seen in the movie of Fig. Efa), is hard to 
occur. In Fig. [2j the enlargement of vortex clusters oc- 
curs in the upper row incidentally. Which of the upper 
and lower rows shows this behavior depends on the small 
noise added to the initial state. 

From the movies of Figs. [T] and [2j we find that the BvK- 
like vortex street emerges faster in Fig. [2] than in Fig. [TJ 
i.e., the vibration of the obstacle seems to enhance the 
emergence of the BvK-like vortex street. In the movies 
of Fig. [TJ vortex dipoles are first generated and change 
into the BvK-like vortex street, which takes some time. 
On the other hand, in the movies of Fig. [2j the BvK- 
like vortex street is generated immediately after the first 
disturbance. 

The behaviors in Fig. [2] is due to the resonance be- 



tween the vibration frequency and the vortex shedding 
frequency. For other vibration frequencies, such behav- 
iors cannot be seen and the BvK-like vortex street is only 
disturbed randomly. 

IV. CONCLUSIONS 

We have investigated dynamics of quantized vortex 
shedding from an obstacle moving in a superfluid, in 
which the obstacle traces a sinusoidal trajectory. As in 
the case of the straight motion of an obstacle, quantized 
vortex-antivortex pairs are released when the velocity 
exceeds some critical velocity, and for appropriate pa- 
rameters, a BvK-like vortex street emerges in the wake. 
In contrast to the case of the straight motion, the sym- 
metry between the two rows of the BvK street is broken; 
the size of the corotating vortex clusters is enlarged in ei- 
ther row. This phenomenon occurs when the oscillation 
frequency of the obstacle is close to the vortex shedding 
frequency. 

The BvK-like vortex street generation in superfluids is 
an example of a highly nonlinear phenomenon far from 
the equilibrium. When some external vibration is applied 
to such a system and the frequency approaches some res- 
onance, we expect that the system exhibits nontrivial 
resonant dynamics due to nonlinearity and nonequilib- 
rium. 
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